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Abstract. Using a subset of the current Run II data, the CDF and D0 have performed several measurements on heavy flavor 
production. In this paper, we present a new measurement of prompt charm meson production by CDF. We also report the latest 
CDF II measurements of inclusive J '/*F production and fo-production without requirement of minimum transverse momentum 
on J/^if and b-quark. They are the first measurements of the total inclusive and b quark cross section in the central 
rapidity region at a hadron collider. The results of J p¥ cross section as a function of rapidity, and fc-jet production cross 
section measured by D0 are also reviewed. 



INTRODUCTION 

Measurements of the production cross section of heavy 
flavor quarks (c and b quarks) in the pp collisions pro- 
vide us an opportunity to test the prediction based on 
the Quantum Chromodynamics (QCD). No only is QCD 
one of the four fundamental forces of the nature, but also 
many searches of new physics beyond Standard Model 
require a good understanding of the QCD background. 

The Bottom quark production cross section has been 
measured by both the CDF and D0 experiments 1 1 , 
2] at the Fermilab Tevatron in pp collision at i/s = 
1.8TeV, and found to be about three times larger than 
the next-to-leading order (NLO) QCD calculation |3, 
4]. Since then, several theoretical explanation has been 
proposed to solve the disagreement, such as the large 
contributions from the NNLO QCD processes, possible 
contribution from "new physics" |5] etc. Recently, a 
more accurate description of b quark fragmentation has 
reduced the discrepancy to a factor 1.7 |g]. Nevertheless, 
further experimental measurements on the heavy flavor 
production are essential to shed light on the remaining 
clouds of this long standing issue. 

Since the end of Run I (1996), the Tevatron had un- 
derdone major upgrade. The center of mass energy in the 
pp collision was increased from yfs = 1.8TeV to y/s — 
1.96TeV. The large luminosity enhancement dramati- 
cally increases the discovery reach and moves the exper- 
imental program into a regime of precision hadron col- 
lider physics. At the same time, both CDF [7] and DO 1 8] 
experimental collaborations significantly upgraded their 
detectors to enrich the physics capabilities, especially for 
the heavy flavor physics program. The Tevatron Run II 
data taking period started in April 2002, so far the ac- 



celerator has delivered about 300pb integrated lumi- 
nosities, both experiments collected more than 200 pb~' 
physics quality data. In this paper, we will give a brief 
summary of the latest results on heavy flavor production 
from Tevatron Run II. 



PROMPT CHARM MESON 
PRODUCTION CROSS SECTION 

In order to help understanding the discrepancy between 
measured Run I Bottom quark cross section measure- 
ment and theoretical calculation, one can repeat previ- 
ous measurement with better control of experimental un- 
certainties. An alternative approach is to examine the 
production rate of another heavy quark. Charm meson 
production cross sections have not been measured in pp 
collisions and may help to solve this disagreement. For 
CDF II upgrade, one of the crucial improvement with re- 
spect to the CDF I is the implementation of a Silicon Ver- 
tex Trigger (SVT) |9], which allows us to trigger on dis- 
placed tracks decayed from long-lived charm and bottom 
hadrons. With this trigger, large amount of fully recon- 
structed charm meson has been collected, opening a new 
window for heavy flavor production studies at hadron 
collider. 

Using just 5.8 ± 0.3 pb~' data, CDF performs a mea- 
surement of prompt charm meson production cross 
section The charm mesons are reconstructed in 

the following decay modes: D° — > K~n + , D* + — ► D°7t + 
with D Q ^R-n + , D+ -^K-n + n + , D+ -> <j>n + with 
<p -^K + K~, and their charge conjugate, as shown in 
Fig.^ We separate charm directly produced in pp in- 
teraction (prompt charm) from charm decayed from B 
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FIGURE 1 

invariant mass distribution of D° 
mass difference distribution of D* 
invariant mass distribution of D + - 
invariant mass distribution of D + 
candidates. 



Charm signals summed over all pj bins: (a) 
-» K~ K + candidates; (b) 
— * D°7t + candidates; (c) 
K~tz + tz + candidates; (d) 
-> (j)7t + and £>+ — > 



mesons (secondary charm) using the impact parameter of 
the charm candidate with respect to the primary interac- 
tion point. Prompt charm meson points to the beamline 
and its impact parameter is zero. The secondary charm 
meson candidate does not necessarily point back to the 
primary vertex, it has a rather wide impact parameter 
distribution. The shape of the impact parameter distri- 
bution of secondary charm is obtained from a Monte 
Carlo simulation. The detector impact parameter resolu- 
tion is modeled from a sample of — > n + n~ decays 
that satisfy the same trigger requirement. The prompt 
charm fraction is measured as a function of charm me- 
son pr- Average over all pj bins, (86.6 ±0.4)% of the 
D° mesons, (88.1 ± 1.1)% of D*+ (89.1 ±0.4)% of D+, 
and (77.3 ±3.8)% of Df are promptly produced (statis- 
tical uncertainties only). 

The measured prompt charm meson integrated cross 
section are found to be a(D°,pT > 5.5GeV/c, \y\ < 
1) = 13.3 ±0.2± 1.5Mb, a(D* + ,p T > 6.0GeV/c, |y| < 
1) = 5.2 ±0.1 ±0.8 Jib, a(D+,p T > 6.0GeV/c,|y| < 
1) =4.3±0.1 ±0.7 jub and a(D+,p T > 8.0GeV/c, |y| < 
1) = 0.75 ± 0.05 ± 0.22 /ib, where the first uncertainty 
is statistical and the second systematic. The measured 
differential cross sections are compared to two recent 
calculations II ll fl 211 . as shown in Fig. |2] and Fig. [3] 
They are higher than the theoretical predictions by about 
100% at low pj and 50% at high pj. However, they are 
compatible within uncertainties. The same models also 
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FIGURE 2. The measured differential cross section mea- 
surements for \y\ < 1, shown by the points. The inner bars 
represent the statistical uncertainties; the outer bars are the 
quadratic sums of the statistical and systematic uncertainties. 
The solid curves are the theoretical predictions from Cacciari 
and Nason ITlll . with the uncertainties indicated by the shaded 
bands. The dashed curve shown with the D* + cross section 
is the theoretical prediction from Kniehl [ 12]; the dotted lines 
indicate the uncertainty. No prediction is available yet for Df 
production. 
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INCLUSIVE J/y PRODUCTION CROSS 
SECTION 

The CDF II detector has improved the dimuon trigger 
system with a lower threshold of 1.5GeV/c. This allows 
us to trigger the J/*¥ — » jl~ H + event with transverse mo- 
mentum all the way down to pj(J /*¥) = OGeV/c. Using 
about 39.7 pb" 1 data, CDF reconstructed 299800 ±800 
J/*¥ candidates (statistical uncertainty only), as shown 
in Fig. A new measurement of the total inclusive Jf¥ 
cross section in the central rapidity region \y\ < 0.6 has 
been carried out. This is the first measurement of the to- 
tal inclusive J/*¥ cross section in the central rapidity re- 
gion at a hadron collider. The differential cross section is 
shown in Fig. [5] and the integrated cross section is mea- 
sured to be: 

a(pp -> //"FX, W/Vl < 0.6) x Br(J/*i> -> /xju) 
= 240 ±l(st at) 



^(sysf)nb, 
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FIGURE 3. Ratio of the measured cross sections to the the- 
oretical calculation from Cacciari and Nason. The inner error 
bar represents the statistical uncertainty, the outer error bar the 
quadratic sum of the statistical and systematic uncertainty. The 
hatched band represents the uncertainty from varying the renor- 
malization and factorization scale. 
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CDF Run II Preliminary 



299,800 ± 800 Events 
Luminosity = 39.7 pb" 1 
Width = 24 MeV/c 2 . 
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FIGURE 4. The invariant mass distribution of triggered J/^f 
events reconstructed in the 39.7pb~' of CDFII data. 



Taking advantage of the large azimuth coverage of its 
muon system, the D0 measured the differential cross 
section of the inclusive J /*¥ as a function of the rapidity 
using 4.74pb~ 1 . The distribution is shown in Fig. [5] 



BOTTOM QUARK PRODUCTION 
CROSS SECTION 

The Run I central b-quark production cross section mea- 
sured by CDF and D0 has a minimum transverse mo- 
mentum cut on the fe-hadrons due to the trigger require- 
ments. Since those measurements only explore a small 
fraction (~ 10%) of the fo-hadron pj spectrum, it is not 
clear from the data whether the excess over the theory is 
due to an overall discrepancy of the bb production rate, or 
it is caused by a shift in the spectrum toward higher pj- 
An inclusive measurement of bottom quark production 
over all transverse momentum can certainly help resolve 
this ambiguity. 

Notice that large fraction of Z?-hadron Hj, decays to 
J/*¥ final states, where H/, denotes both hadron and anti- 
hadron. CDF performed an analysis to extract the in- 
clusive fe-hadron cross section from the measured in- 
clusive J production cross section. Due to the long 
live lifetime of the fe-hadrons, the vertex of the J '/I* de- 
cayed from a fo-hadron (secondary // V F) is usually hun- 
dreds microns away from the primary vertex, while for 
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FIGURE 5. The differential cross section for inclusive J [*¥ 
as a function of pj with |y| < 0.6. 

the prompt //*P's, which are directly produced or de- 
cayed from higher charmonium states, their vertex are 
at the proton-antiproton interaction point. As the result, 
we can statistically separate these two components by 
examining the projected J /*¥ decay distance along its 
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FIGURE 6. The differential cross section for Jf¥ as a func- 
tion of |y|. 
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FIGURE 7. fo-hadron differential cross section as a function 



transverse momentum. The measured ^-fraction is then 
applied to the previously measured inclusive 7/ v P cross 
section to obtain the differential production cross section 
of H b — > J /*¥X as a function of pj(J However, the 
above algorithm to extract ^-fraction does not work for 
y/ 1 ? candidate with low transverse momentum, because 
a fo-hadron with low p? does not travel far away from the 
primary vertex in the transverse plane. Therefore a mini- 
mum pr (J f¥) > 1.25 Ge V jc requirement is imposed for 
the fo-hadron cross section measurement as a function of 
/?r(.// v F) in order to have a reliable determination of b- 
fraction. The measured differential cross section result is 
shown in Fig.Q 
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FIGURE 8. £>-hadron differential cross section as a function 
of p T (H b ) with |y| <0.6. 

Luckily, the fe-hadron decayed at rest can still trans- 
fer a few (~ 1.7) GeV/c transverse momentum to its 
J/*¥ daughter because of the large b-hadron mass. The 
knowledge of the secondary J /\P production with trans- 
verse momenta less than 2.0GeV/c allows us to probe 
the b-hadrons with low transverse momenta down to 
zero. Therefore using the measured fc-hadron differen- 
tial cross section with 1.25 < p T {J < 17.0GeV/c, 
we are able to extract the fo-hadron differential cross sec- 
tion as afunctionof pt(Hi,) down to pr(Hb) — QGe"V/c. 
To do so, we perform a convolution that is based on the 
Monte Carlo template of the fo-hadron transverse mo- 
mentum distribution. The unfolding procedure is then 
repeated using the measured fo-hadron production spec- 
trum as the input spectrum for the next iteration, until 
the x 2 comparison between the input and output spec- 
trum calculated after each iteration becomes stable. The 
measured fr-hadron cross section as a function of pj (Hi,) 
is plotted in Fig. [8] The total fe-hadron cross section is 
found to be 

a(pp -> H h X, \y(H b )\ < 0.6) x Br(H b -> JfifX) 
xBr(J ^ nn) = 24.5±0.5(stat)tf 7 g (syst)nb, 

(2) 

The total single Z?-quark cross section is obtained by di- 
vided this measurement by two, the branching fractions, 
and the rapidity correction factors obtained from MC. We 
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FIGURE 9. Measured Run II £>-jet cross section compared 
to theoretical prediction. 



find: 

a(pp^bX,\y(b)\< 0.6) 

(3) 

= 18.0±0 A(stat)±3.8{syst) fib, 

and 

a{pp^bX,\y{b)\<\) 

(4) 

= 29.4 ± 0.6(stat) ± 6.2(syst ) jib. 
B- JET PRODUCTION CROSS SECTION 

Using 3.4pb~' of data, a b-jet production cross section 
analysis has been performed by D0 . The candidate event 
is selected by associating a muon track with a jet, which 
is defined with R = \J A77 2 +A0 2 cone algorithm. Notice 
the fact that the muon decayed from b quark has a higher 
transverse momentum with respect to the net momentum 
vector of combined muon and jet, the corresponding 
distribution is then used to fit the b-jet fraction. The 
signal template is modeled from a b — > jx Monte Carlo 
simulation, and the background template is extracted 
from 1.5 million QCD events. After correcting the muon 
and jet reconstruction efficiency, and the calorimeter jet 
energy resolution, the fo-jet cross section is obtained. The 
preliminary result is shown in Fig. [9] compared to the 
theoretical prediction. This measurement is consistent 
with the previous Run I measurement 11411 . 



CONCLUSION 

The understanding of the heavy flavor production is cur- 
rently one of the most important challenges faced by 
QCD. In this paper, we present the most recent mea- 
surements from CDF and D0 experiments. Given large 
amount of data that have already been collected, rapid 
improvement of Tevatron performance, and further un- 
derstanding of the detector, we expect that much more 



precise measurements on the heavy flavor production 
will be available in the near future. At the same time, 
several other analyses have also been carried out to ex- 
plore other aspects of heavy flavor production mecha- 
nism, such as the bb correlation and cc correlation during 
the production. In the next a few years, a combination of 
better experimental data and improved theory should ad- 
vance our knowledges of the heavy quarks production. 
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